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Abstract 

Corrections to a mean number of droplets appeared in the process 
of nucleation have been analyzed. The two stage model with a fixed 
boundary can not lead to a write result. The multi stage generalization 
of this model also can not give essential changes to the two stage 
model. The role of several first droplets have been investigated and it 
is shown that an account of only first droplet with further appearance 
in frame of the theory based on the averaged characteristics can lead 
to a suitable results. Both decay of metastable phase and smooth 
variations of external conditions have been investigated. 

1 Introduction 

Up to the last years a kinetic descriptions of nucleation processes were based 
on the averaged intensity of droplets appearance, i.e. on the rate of nucle- 
ation. Namely the rate of nucleation is considered as a central characteristic 
of the first order phase transition. But the supercritical droplets appear in 
the system occasionally and this feature has to be taken into account. Since 
[Hj, [Sj, [Z| several attempts to include stochastic effects of droplets appear- 
ance were made. Unfortunately, they could not give an adequate and rather 
precise description of stochastic effects. 

The arguments in justification of kinetics based on the averaged character- 
istics (see j3]) remain valid even after stochastic formulation of the nucleation 
problem has been presented. When the number of droplets in the system is 
very big the result of the theory based on the averaged characteristics (TAG) 
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is precise. Namely the number of droplets is the central characteristic of the 
nucleation process and it is calculated in experiments. When the number of 
droplets appeared in a system isn't so great (in the free molecular regime 
of growth it occurs only due to a volume of a system) one can speak about 
corrections to a mean value of the total number of droplets appeared in the 
system in comparison with result predicted by TAG. This value will be the 
main object of investigation here. 

In diffusion regime of the droplets growth the kinetic description is based 
on special models (see jT]) and there is no need to analyze this regime here. 
So, in this paper only the free molecular regime of droplets growth will be 
considered. 

In investigation of stochastic effects of nucleation one can see rather easy 
that the first correction terms are equal to zero. So, there appear a prob- 
lem to calculate the higher terms of decomposition. It is rather difficult to 
perform precise calculations and we shall be interested at least to get es- 
timates for these coefficients to prove the smallness of the total amount of 
corrections. But even the calculation of the very first terms meets technical 
difficulties (see [7], j^]). One has also to stress that the zero shift found in 
these papers was a natural result of linearization made to overcome technical 
difficulties. Then it can not be considered there as a true physical result, 
but only the consequence of lineariazation. So, it is necessary to propose a 
method to calculate the estimates for coefficients in asymptotic expansions 
due to stochastic corrections of nucleation. 

There are two characteristic situations of external conditions in which 
kinetics of nucleation ordinary was constructed. These conditions are: 

• decay of metatsable phase when at some moment the metastable phase 
is created and later there is no external influence on the system; 

• smooth variation of external influence on the system. 

In both situations corrections will be established. 
The structure of the further analysis is following: 

• At flrst we shall analyze the two stage model with a flxed boundary. 
The result will be disappointing - one can not reproduce the results of 
numerical simulation. This corresponds to the difficulties of this model 
in prediction of the value of dispersion. 
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• To improve results we shall use the multi stage generalization of the 
last model. But corrections to the two stage model will be small and 
this can not lead to suitable results. 

• Then a new approach will be used. We shall analyze the role of stochas- 
tic appearance of the first droplets. Results will be very fruitful and 
one can see that already account of one droplet will lead to success. 

• All these considerations will be made both for decay of metastable 
phase and for the smooth variations of external conditions. 

2 Decay of metastable phase 

The kinetics of nucleation in frames of the theory based on averaged charac- 
teristics (TAG) can be described by the following equation 

g(z) = / (z — x)^ exp(—g(x))dx 
Jo 

where unknown function g is the renormalized value of the number of molecules 
in a liquid (new) phase. This result can be found in jH]. The meaning of vari- 
ables z, X can be also found in |H]. Since jH] it is known that one can describe 
kinetics in frames of monosdisperce apporoximation, i.e. 

9iz) = Neff{z)z' 

where A^e// is the effective number of droplets in monodisperce peak, namely 

N,jf{z) = z/A 

The monodisperce approximation can be chosen as the fixed (not floating) 
monodisperce approximation (see |H1) and leads to the following expession 
for the size spectrum 

f{x) = f,exp{-NeffX^) 
Here /* is the amplitude of spectrum, 

N^ff = NiAx/A) 

and Ax is a width of a whole spectrum (connected with the duration of a 
nucleation period). 
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The total number of droplets can be obtained on the base of / as 

Ntot = / dxf{x) 
Jo 

For Neff we have a Gaussian distribution with standard dispersion since 
formation of the first Neff droplets can be treated as the sequence of inde- 
pendent events 

^ 2 < Neff > ' 

Here < N^ff > is the mean value of N^ff 

Then for the averaged value of Ntot, i-e. for < Ntot > we have the following 
formula 

/oo roo 
dNeffP{Neff)U / e^v{-NeffX^)dx 
-oo Jo 

Now with the help of formula 

dy exp{—ya) exp{—c{y — 6)^) ~ exp(— 6a + — ) 

we fulfill integration over N^ff. Here b —< Neff >, a — x^, — 2 < 
Neff >. As the result we have 

/•oo ^6 

< Ntot >~ / exp(- < Neff >x^ + — < Neff >)dx (1) 
Jo I 

The second term in exponent, i.e. ^ < Neff > is the correction term which 
can be seen from 

< Ntot >~ 0' HeM-y' + ^^t J ^y 

Jo 2 < Neff > 

As a rough estimate we can take y in correction term as y ~ 1 and get 

< Ntot >=< Ntot > exp(-— — ^ -) 

2 < Neff > 

where < Ntot o > is the value Ntot calculated without stochastic effects taken 
into account, i.e. in frames of TAG. Then for this value one can get expression 

< Ntot >~ exp(- < Neff > x^)dx 

4 



Having noticed that 

< Neff >^< Ntot > /4 

we get 

2 

< Ntot >=< Ntoto> exp 



< Neff > 
Decomposition of exponent gives 

2 

< Ntot >=< Ntot > +2 + — + ... 

< Ntot > 

Another more balanced variant of consideration is to use decomposition 
of g 

exp(y < N,jf >) = l + ^<N,ff>+^< N,jf 

already in At least the integral then will have no problems with conver- 
gence. We have 

< Ntot >~ / exp(- < Neff > x^){l + — < N,ff > +— < N^ff >^)dx 
Jo Z o 

Integration can be fulfilled separately for every term in decomposition. Then 
we come to 

< Ntot >=< Ntot > +^1 + — TT + ... 

< Ntot > 

where constants Ai and A2 are given by 

A, = 2 exp{-y^)y'dy/ / exp{-y^)dy = 8/9 
Jo Jo 

roo roo 

A2 = 2 exp{-y^)y^^dy/ / exp{-y^)dy = 510/89 = 6.91 
Jo Jo 

From the functional forms of expressions for Ai one can see that Ai is 

determined rather smart while the error in A2 can be essential. The reason is 

the rapidly growing term y^^ in subintegral function. Already y^ in expression 

for Ai grows too rapidly. So, in subintegral functions the main role belong 

to y corresponding to droplets appeared at the very end of the nucleation 

period. But the form of spectrum ~ exp(— y^) is determined at the back side 

of spectrum (i.e. at y > 1 ) with a low accuracy. In TAG the weight of such 

droplets was negligible and the result was accurate. Here the error can be 
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essential. That's why it is reasonable to restrict the decomposition only by 
the first term Ai. 

In Figure 1 one can see the results of numerical simulation (oscillating 
curve) and analytical approximation (smooth monotonuous curve) for the 
relative value 

P =< Ntot > / <Ntoto> -1 
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Figure 1: Numerical and analytical solutions in the situation of decay. The 
initial monodisperce approximation is considered. 

One can see that there is no satisfactory coincidence between the theo- 
retical result and the result of simulation. The reason is the roughness of 
monodisperce approximation used in [H] and applied here. 

Now we shall take a more refined approximation used to calculate the 
value of dispersion initiated by stochastic appearance of droplets in the pro- 
cess of decay |9i|. This approximation is the following: the length of formation 
of monodisperce spectrum is 2 * / where / = 0.2; the monodisperce spectrum 
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is formed at 2 * Z — 6 where b = 0.336. The derivation of this approximation 
can be found in jOj. Here 

POO 

Jo 

Then one has to reconsider the value of Ai. It can be done only rather 
approximately. One can see that Ai is reciprocal to the total number of 
droplets in monodisperce spectrum which is now 2*1 = 0.4 instead of 0.25 
in initial monodisperce approximation. Then instead of previous Ai one has 

to take 

Ai AiO.25/0.4 = 1/3.6 
The results are shown in Figure 2 
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Figure 2: Numerical and analytical solutions in the situation of decay. The 
shifted monodisperce approximation is considered. 

Here the coincidence between curves became better but it is not still 
satisfactory. This means that the model with a fixed boundary can not give 
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a good result. This corresponds to the difficulties in calculation of dispersion 
on the base of the model with a fixed boundary [TTj . 



3 Smooth variations of external conditions 

In the situation of smooth variation of external conditions (so called dynamic 
conditions) we can fulfill the same procedure. 
The evolution equation in TAG looks like 

g = / (2 — x)^'iIj{x) exp{—g{x))dx 



where = exp(x) describes the change of external conditions and the renor- 
malization to cancel the coefficient in is used. As a compensation for such 
renormalization the coefficient c = 0.189 or c = 1/6 (it depends on the type 
of choice of the base of decompositions (see jl|)) appears. 

Here we shall use the monodisperce approximation also. The monodis- 
perce approximation was proposed in |10| where all details can be found. In 
the fixed monodisperce approximation one can write 

f = f,exp{x-<Neff>{x + 3f) 

for the spectrum of droplets sizes. Here the effective number of droplets is 
given by 

Neff = ^, c= 0.189 

or 

Neff = — 

en 27 

For the total number of droplets we have an evident expression 

/oo 
dxexp{x- <Neff>ix + 3)3) 
-00 

For the mean total number of droplets one can get 

//■oo 
dNeff / dxexp{x- < Neff >{X + 3f)P{Neff) 
J — 00 

where P^N^ff) is the distribution function for the quantity of effective droplets. 
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For P{Neff) we have an evident Gaussian distribution with dispersion of 
ideal gas 

P(Neff) ~ exp(- ^ " 

Z < I\eff > 

Having fulfilled integration one gets 

/oo f/j. _|_ 3~)6 

dxexp{x- < Neff > {x + Sf) exp( ^' < TVe// >) 

-oo Z 

Certainly, the last integral doesn't converge. We need a regularization 
which will be done below. 

Then one has to decompose 
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exp( ^ ^ < Neff >) = 1+ ^ ^ < Neff > g < Neff >' +... 

and then one can fulfill integration for every term. Now the integral has no 
problems with convergence. 

In above formulas < N^ff >= 1/27 which is a natural requirement to use 
the monodisperce approximation. To calculate correction terms one has to 
include explicitly the volume of the system V (i.e. the real mean number of 
droplets). Now we shall give the exphcit formulas for correction terms. We 
have to calculate the value 



< N,^, >= r dx r exp(x -^{x + 3f) eM-^-^^^^J^dNeff 
J—oo J—oo V ZrJgff 



Neff =< Neff > 

Earlier we decomposed exp(— ^^''^^~^^''J^^ ) and had some problems with con- 

V 



vergence. Now we shall decompose exp(x v^{x + 3)^). At first we shall 



present this exponent as 

exp(x - ^{x + Sf) = exp(x - ^{x + Sf) exp(-(^ - ^){x + 3f) 
The decomposition of the last exponent gives 

, .Neff Neff., ^ , ^v^Neff - N~eff (x + 3f (Neff - N^fff 

exp(-(-^ ^){x+3f) = l+{x+3f ^ '^^ + ^ ^ ' ^ + 
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6 \/3 24 

The calculation of integrals gives 

We have to notice that N^ff — V/21. Then we have the decomposition 

, , 1 ~ r^eic£>ix- ^{x + ?>f){x + ?>fdx 

exp{x -^{x + 3y)dx 

1 ~ , f^^expix-^{x + 3r){x + 3y'dx 

Now we have to note that the lower limit of integrations has to be put 
X = —3 because the monodisperce approximation begins to work only at 
X > —3. The region a; < 3 has negligible influence in the total amount of 
droplets. Then 



< Ntot >=< NtotiV = oo) > (1+ 



1 2NlffJ^3exp{x-^{x + 3f){x + 3fdx 
J^^expix- ^{x + 3f)dx 



3 ~ 2 l^,exp{x-l^{x + 3f){x + 3y^dx 
+ 241/4 ^ff jo^exp(x- ^(a: + 3)3)dx 

Here there were no problems with convergence. We can calculate the 
integrals numerically which gives 

< Ntot >=< NtotiV = oo) > {1+ —+ 

< NtotiV = oo) > < NtotiV = oo) >^ 

1 JZoCxpix-^ix + 3f)ix + 3)^dx _ 
' 2*272 J^^expix-^ix + 3y)dx 
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' 8*274 J^^exp{x-^{x + 3y)dx 

This is the final result. 

Again one has to note that A2 is determined with uncertainty caused by 
approximate description of the back side of spectrum. Here A2 is calculated 
to see that there is no singularities in decomposition. Again it is reasonable 
to take into account only the first correction term. 

In the situation of the smooth behaviour of external conditions there is 
no moment of start. The point of formation of monodisperce spectrum is 
z = —3 and it is determined in the internal point. So, when we observe 
the subintegral function it has a maximum not in the boundary point (as in 
decay when it is 2; = 0) but in an internal point near z = —3. The subintegral 
function {z — x)^il!{x) exp{—g{x)) is rather symmetric around x = —3 (in 
decay one can not imagine that the subintegral function is symmetric around 
the boundary point). 

Numerical simulation and analytical result can be seen in Figure 3. 

One can not see the satisfactory coincidence between theoretical result 
(monotonuous curve) and numerical siulation (oscillating curve). So, the 
used model can not give the good results. 

To the data of numerical simulation one can suggest a phenomenological 
approximation 

Ar~F(l + ^ + ...) 
A = 5/6 

One can establish the universality also in stochastic formulation of the 
problem. Then the parameter A is no more than a universal constant. The 
numerical simulation gives A = 5/6. Generally speaking one can stop here 
all investigations. 

So, the model with a fixed boundary failed in determination of correc- 
tions to the number of droplets. One has to turn to the models with fioating 
boundary which was succesfuUy applied in jU], [T2] to the calculation of dis- 
persion. But in the models with a fioating boundary the correction to the 
mean number of droplets is zero. This isn't the error of the model but the 
level of description is limited here. 

The next step in consideration is to analyze whether one can reconsider 
the previous results . The fact is that in nucleation one can see some balanc- 
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Figure 3: Numerical simulation and analytical solution under the smooth 
external conditions. The dependence of P over the volume of the system 
V = X is plotted. 

ing forces. This can require the addition more detailed analysis containing 
several stages model instead of two stages. 

4 Three stage scheme 

In nucleation under the smooth behavior of external conditions one can ob- 
serve the specific property of compensation which can damage the previous 
consideration. One can analytically observe the fact of compensation briefly 
described below which means that one can not use two stage scheme to get 
corrections for the mean number of droplets. The result of the three cycle 
scheme can be the leading term in the shift of the mean droplet number. So, 
we have to use the three stage scheme. 
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The effect of compensation in the two stage scheme require to give esti- 
mates in a three stage scheme. At first we shall explain the cancellation of 
eff'ects in the two stage scheme. 

The number of droplets Ntot can be in TAG in the monodisperce approx- 
imation calculated as 



/oo 
-3 



Here N^ff is 1/27. 

If due to stochastic effects the necessary number N^ff appeared up to the 
"moment" z — —3 — 5 instead of 2; = — 3 then we have to calculate Ntot as 

/oo 
exY>{x- N^ff{x + ?, + 5f)dx 
-3— <5 

The calculation gives 

/oo 
exp(-5) exp(|/ - N^^^^y + 3 + 5f )dx 
-3 

ior y — X -\- 5 and 

Ntot{5) =< Ntot > exp(-(^) 

Now we shah estabhsh the distribution P{5) of the shift 5. The distribu- 
tion P{Ngff) is the ordinary Gausiian distribution 

To get P{6) we use 

P{N,ff)dN,ff = P{S)dS 
The derivative dN^ff/dS is 



dNeff 



exp((5) 



d5 
Then 

P{5) = P{N,ff) exp(5) 

The second factor completely compensates the shift in the total number of 
droplets. Really, 



< Ntot >= J P{S)Ntot{S)dS = TV, 



mean 
tot 
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This compensation shows the zero effect in the shift of droplets number 
in the two cycle scheme and requires to consider the three stage scheme. 

Consider the three stage scheme qualitatively. We have to mention that 
the value A^e// doesn't purely appear under the ideal conditions. Already at 
z — —3 the small part of substance is in the droplets. The main consumers 
of vapor at ^ = — 3 are the droplets appeared at z — —6. Here we can use 
also the modisperce approximation and present g at z — —3 as 

g{z = -3) = Nir^itiz + 6)^ 

with parameter Ninu of initial monodisperce approximation. 

Since the stochastic number Ninu doesn't coincide with the value Ninu 
calculated in TAG we can see the deviation of stochastic value g{z — —3) 
from the value g{z = —3) calculated in TAG. 

It seems that we come to the situation which has been already described 
in the two stage scheme. But now we don't observe the effect of compensation 
because this effect takes place only due to the integration in the infinite limits. 
But here such an integration is absent - moreover we need the effect up to 
the fixed moment z = —3. So, we shall see the effect which results in the 
difference of the mean value < g{z = —3) > from g{z = —3). Then we see 
the regular shift 6z of the moment until which the number of droplets A^o 
appears. This regular shift leads to the regular shift in the total number of 
droplets. 

Now we fulfill the computations. The number of droplets formed until 
zi = —3 in TAG is 

/■- 3 

^eff = / exp(x - kNinit{x + 6f)dx 

J —oo 

Parameter k has here a role like in the two cycle scheme had. 

The distribution PiyNinit) of the stochastic number of initial droplets Ninit 
has a normal Gaussian form 

P(7V,„,,) ~ exp(- ^^--:~^--^^' ) 

The value of dispersion here corresponds to the fact that we have a free 
stochastic appearing of droplets. 

Here appeared a special question whether it is possible to write the gaus- 
sian distribution for the total number of events (here it is the number of 
appeared droplets). But this question can solved positively. 
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Then 

< Neffiz = -3) >= J P{Ninit) exp(a; - kNinit{x + Qf)dx 

Hence 

< ^e//(^ = -3) >^ N^ff 

Then we can calculate the regular shift Sz along z-axis for the moment 
until Neff droplets appeared. For the value of 5z we get 

exp(— 3)- 



5z- — ^ 



exp(x|a;=-3 - kNinit{x\a;=-3 + 6)3) 

Denote by S the following value 
Let US calculate < Neff > 

< N,ff{z = -3) >= r dS r' dx^^^^^^£^^^^eM^-k{Ninit+5){x+6f)V 

Recall that V is the volume of the system and here we have to introduce it 
explicitly. 

We fulfill the calculations and get 

< Kff{z = -3) >= 
d5 r'dx^^^^^^£S^^£^eM^-kN,r,it{x+6f)V{l-^^^ 

Then the deviation between < Neff{z — —3) > and < Neff{z — —3) > 
will be 

< Neff{z = -3) > - < Neffiz = -3) 

1 roo r — 3 _ _ 

dy dxeM-y^) exp(x - kNi^uix + &T)—y\x + 6)%2Ni^uV)-^ 

\/7r J —oo J —oo Z 



Since the monodisperce approximation becomes suitable only at 2; = — 6 
it is more reasonable to write 

< Neffiz = -3) > - < Neffiz = -3) >= 
15 



— / dy dxexp{-y^) exp(x - kNi^u{x + Qf)—y\x + 6f2N,nitV-^ 

yTT J-oo J-6 Z 

Then for the average number of droplets appeared in the system 

< N,ff{z = -3) > y = {Neffiz = -3)U^^^(,=_3)=exp(-3) + CoV-')V 

where 

][ roo /•— 3 - A;^ - 

Co = ^ / c?2/ / (ia;exp(-|/2) exp(a; - A;iVi„it(a; + Q)^)—y^{x + 6)^2iVi„it 

yTT J-oo J-6 I 

Here we use N^ff^z = —3) = exp(— 3) for the number of droplets in the unit 
of the system because namely this value corresponds to the characteristic 
unperturbed value of droplets and the choice 

Nir,it = exp(-6) 

which was used in the last formula. 

The numerical calculations gives the value of shift Sz 

Co 0.01 

dz — ~ 

exp{-3)V V 

The shift in the total number of droplets can be easily calculated since 

< Ntot >= exp{Sz)Ntot 

Hence 

- w 
w ~ 0.01 

Quite analogously one can show the smallness of corrections in a three 
cycle scheme in the situation of decay. 

The result of performed calculations shows the smallness of corrections 
to the mean value of the total number of droplets appeared in the processes 
of nucleation under the conditions of decay and under the smooth behavior 
of external conditions. 

Correction terms calculated by the theoretical derivations in the three 
cycle scheme are so small that they can not be confirmed both by numerical 
simulations and experimental researches. So, generally speaking one has to 
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say that these stochastic corrections don't appear in any practically signifi- 
cant terms of decompositions. 

This result shows that there is no sense to fulfill the procedure of renor- 
malization analogous to that used in ^1] in calculation of dispersion in terms 
of two cycle explicit model with a fixed boundary. 

Estimating the total result we have to stress that the weak feature of the 
presented method is the approximate knowledge of the back side of the size 
spectrum exp(— x^) in decay and exp(a; — (x + 3)^/3^) in the smooth variation 
external conditions. So, these constructions can not be appreciated as the 
concretely determined result. 

5 Several first droplets 

The leading correction term doesn't depend on the volume of the system. 
It means that the deviation in the mean number of droplets from the value 
predicted by TAG doesn't increase with increase of the volume. The latter 
means that namely the first several droplets are responsible for initiation of 
corrections to the mean value of droplets. So, now we shall show how many 
droplets are responsible for these corrections. 

We know that on one hand the Gaussian distribution can not be applied 
as statistics for the first droplets and on the other hand the big droplets are 
extremely important in kinetics, which lies in the base of iteration method. 
So, it is reasonable to use the explicit numerical simulation to see the role of 
several first droplets. 

We start consideration of the role of several first droplets with the case 
of decay. Figure 4 illustrates the role of the stochastic appearance of the 
first droplet in nucleation kinetics. Here the relative excess of mean droplets 
number is shown. 

There are two curves, both are functions of the volume of system V. 
The value of V is connected with the total number of droplets in TAG as 
Ntac = 1.28V. 

The broken line is the result of numerical simulation for initial problem, 
the smooth line is the result of solution of the following problem: The first 
droplet appears stochastically and later all other droplets appear with prob- 
ability 

pdx ~ I{—)dx 
ax 
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Figure 4: Relative excess of the mean droplet number. Situation of decay 

Here / is the rate of nucleation. So, except the first droplet the further 
appearance occurs according to TAG. 

One can see the satisfactory coincidence between the model and the sim- 
ulation of initial problem. 

One has to note that we are interested in corrections to the droplets num- 
ber when they are essential. We are not interested in the tails of asymptotics. 

The next picture illustrates the model with two stochastically appeared 
droplets. The broken line is the numerical simulation and the smooth line is 
the model with two stochastically appeared droplets. 

Here the coincidence between the model and simulation is practically per- 
fect. But the model with the first stochastically appeared droplet is suitable 
also and due to simplicity has to be considered as the basic theoretical model 
explaining the corrections to the mean number of droplets. 

One can also investigate the model with discrete regular appearance of 
droplets. One can adopt that all droplets appear when 



attain integer values. In this model one can take that the first droplet appears 
stochastically. Nothing will be changed. 
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Figure 5: Relative excess of the mean droplet number. Situation of decay. 
The model with two stochastically appeared droplets. 

The results are shown in figure 6. The axis are the same. One can see 
that the lower broken line which is the result of the last model has nothing 
in common with the upper line which is the result of simulation. So, we has 
to conclude that the discrete effects don't manifest themselves in nucleation 
kinetics. 

One can observe one interesting feature of kinetics. If in the first moments 
of nucleation period the number of appeared droplets is higher than the 
average value then the total number of droplets will be lower than the average 
total value of the droplets number. This effect will take place at rather big 
value of the total number of droplets. At the small numbers of the average 
total number of droplets the effect will be the opposite one. 

Now we shall turn to investigation of the nucleation under the smooth 
external conditions. 

Figure 7 shows results of regular continuous solutions with several first 
droplets born stochastically. There are three curves drawn in this figure. The 
oscillating curve is numerical solution, two smooth curves are approximations 
with the only first droplet born stochastically and with the first two droplets 
born stochastically. It is clear that there is no big difference between these 
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Figure 6: Relative excess of the mean droplet number. Situation of decay. 
The discrete model and simulation of initial problem. 

curves. It means that it is sufficient to take into account only the stochastic 
appearance of the first droplet. 

We see that the coincidence is rather satisfactory. 

Now we shall see how one can incorporate the discrete effects in this situ- 
ations. We propose the following model. The first droplet appears stochasti- 
cally and later droplets can appear only after the elementary fixed intervals. 
Every interval is chosen to have the integral of the rate of nucleation over 
time equal to V~^. The vapor is consumed by the finite (big) number of 
droplets born in the mentioned moments of time and growing regularly. The 
result at small V is shown in Figure 8. 

We see that the coincidence is satisfactory. May be it is even better 
than the result of the regular continuous model with a first stochastically 
appeared droplet. In any case we see that the stochastic appearance of the 
first stochastically droplet diminishes the role of discrete effects. 

One can see that the deviation between discrete and continuous models 
is not big, moreover we see that the role of discrete effects is not essential in 
the deviation of the average number of droplets. 

In this point the nucleation under conditions of decay differs from the 
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Figure 7: Numerical solution and stochastic approximations at small V. 
Smooth external conditions. 

nucleation under the smooth behavior of external conditions. 

The main result of performed simulations is that the stochastic deviation 
of mean value of droplets is mainly caused by the stochastic appearance of 
the first droplet. The stochastic appearance of the first droplet is very simple 
to calculate analytically. Really, we have to write the Poisson distribution 
for the probability of appearance of the first droplet 

P ~ exp(— /) 

or more concretely for staying without appearance of any droplet. Here / is 
the number of possible events. Now we have to come from / to the time t. 
This connection is given by 

/ = exp(t) 

in appropriate renormalization of time t. Certainly this connection corre- 
sponds to the linearization of ideal supersaturation (see |4 ) and the ideal rate 
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Figure 8: Numerical solution and stochastic discrete approximation at small 
V. 

of nucleation will be like exp(x). Until the appearance of the first droplet 
the supersaturation is certainly the ideal one. Then it is easy to get the 
differential distribution p over t as 

pdt = Pdl 

Then 

p = exp(t) exp(— /) = exp(t — exp(t)) 

It is remarkable that the last distribution is the same as the universal 
distribution of droplets in TAG established in [Hj. 

One can write the Poisson distribution for the appearance of the first 
droplet pi ~ exp(— /), for appearance of the first two droplets p2 ~ exp(— Z)/2, 
for appearance of n droplets p„ ~ exp(— /)/n!, etc. When n <^ Ntot one 
can use / = exp(t) to recalculate p{t) on the base of P{1)- This restriction 
isn't essential because at least N^jf <^ Ntot and one needs n < N^jf. 



22 



The effects of discrete model can be also described analytically. It is 
simply necessary to substitute in TAG the integral by the sum. One can act 
in two ways. 

The first possibility is to take explicitly into account several first droplets 
(let it be K). Then the number of droplets in a liquid phase will be 

9 = ~ ^i)^ + f* - exp(x - g{x))dz 
i -^^c 

where is the "amplitude of spectrum" (see [2]) and Zi = Xi are determined 

by 

/* / exp{x)dx = exp(zj) = {i + 0.5) 

J —oo 

and 

I exp(x)(ix = exp(2;c) = (^ + 1) 



Then the methods of solution are quite analogous to 

Another possibility is to use the Euler-McLorrain decomposition for 



zZ^i ~ exp(a;)(ia: 



This approach leads to Bernoulli numbers and will be published separately. 

Also one can use use discrete approximation for all droplets and replace 
it by the integral with the help of the global Euler-McLorrain decomposition. 



6 Concluding remarks 

Generally speaking the most important result of the given consideration is 
the zero correction in the main term of the shift of the mean droplets num- 
ber. The consequence is the conclusion that only several few droplets are 
responsible for corrections in the mean number of droplets. The fact that 
only several first droplets form correction in the total number of droplets 
is rather important for applicability of proposed method to calculate cor- 
rections. The use of monodisperce approximation with a fixed boundary is 
possible only in the case when the first correction term is the zero one and 
the first nonzero term corresponds to the finite (independent on V) absolute 
shift in the number of droplets. Only then the shift is initiated by several 
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first droplets and there is no difference whether we take into account the rest 
droplets in the monodispcrcc peak or not. One can note that the difference 
between the fixed boundary and the fioating boundary is reduced only into 
account of the rest droplets. So, there is no difference what type of boundary 
is used (this isn't true for other characteristics like dispersion). 

Nevertheless we shall give the corresponding derivation in frames of float- 
ing boundary. This is done to show the role of non-gaussian effects (the dis- 
tribution isn't the gaussian one). Having written the expression for < Nfot > 



J^3+y dxex.p{x - ^{x - (3 + y))^)dx 



where y is the shift and P{y) is the partial distribution over coordinate y, 
one can get corrections for < Ntot >■ 

For the partial distribution p{y) one can write 

p(y)dy = P(N)dN 

where P{N) is the partial distribution over possible droplets A'^. Having 
written 

dN/dy = dex.p(N)/dy = exp(A^) 
one can get P{N). For P{N) one can write the ordinary gaussian distribution 

(N- < N >Y 
P{N) ~ exp(- ^ ' ' ^ ) 

After transformations we see that 

(exp(-3 -y)- exp(-3))\ 



P{N) = exp( 



2 exp(-3) 



or 

^ ' 2exp(-3) ' 

Having fulfill decompositions we get 

^ exp(-6)(l-, + ,V2 + ...-l)^ 

^ ^ 2exp(-3) ^ 

and with restriction of first terms 
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or finally 




exp(— 6)?/^ 
2exp(-3) 



Having extracted Gaussian distribution we get 



P{N) = exp( 



exp(— 6)?/^ 
2exp(-3) 



)(1 



exp(— 6)7/^ 
2exp(-3) 



+ ...) 



So, there appear the non-gaussian corrections. Namely these corrections are 



the reason of appearance of corrections in the total number of droplets. They 
can be easily calculated by the manner described in calculations in the model 
with a fixed boundary. 

Here we shall stop our calculations and put a question what distribution 
has to be a gaussian one: the distribution P{N) or the distribution P{y)'? 
Certainly there is no clear answer on this question. Moreover the results of 
[TT] shows that there is a real difference when gaussian distribution instead 
of the Poisson distribution is used. Here it is clear that we have to use 
the Poisson distribution. But then to fulfill the integration one has to use 
the steepens descent method which is equivalent to the use of the gaussian 
distribution with corresponding corrections. 

So we came to a paradox and it can not be resolved without taking into 
account that several first droplets are the reason of the shift of the mean 
number of droplets. Fortunately, there is no need to continue this procedure 
and one can take into account the influence of several first droplets explicitly 
by the procedure described above. 

In investigation of the shift to the droplets number one has to take into 
account that the asymptotic we need is the "intermediate asymptotic" . There 
is no necessity to know for example that instead of 10000 there will be 10005 
droplets. We need the shift where it is at least few percent. So, we need 
asymptotics at the intermediate mean number of droplets. Namely this case 
was investigated and it was shown that already account of two or three first 
droplets is sufficient for the true shift of the droplets number. 

We have to note that there is another reason of applicability of monodis- 
perce approximation with a fixed boundary. This reason lies in construction 
of monodisperce approximation and it is different for decay and for smooth 
variation of external conditions. For the situation with the smooth variation 
of external conditions one can note that the amplitude exp(— 3) correspond- 
ing to the moment of formation of monodisperce peak is very small. So, the 
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value |/exp(— 3) will be small. Then we can neglect y in the lower boundary 
of integration and come to the approximation with a fixed boundary instead 
of approximation with a fioating boundary. So, the smallness of exp(— 3) is 
the reason why one can use approximation with a fixed boundary. 

In the situation of decay there is no such smallness of amplitude. But one 
can recall that the shift of monodisperce approximation (i.e. the position of 
peak formation) was chosen in such a way that "the length" of peak corre- 
sponds to the extremum of droplets number (see [HI). So, the derivative of the 
total number of droplets over the length of spectrum is zero and there is no 
difference whether to use the fixed boundary or to use the fioating boundary. 
Here appears the physical reason of the choice of monodisperce approxima- 
tion in a way prescribed in [D]. Certainly neither the approximation of fixed 
boundary nor approximation of fioating boundary refiect the right physical 
evolution but such a choice of monodisperce approximation allows to ignore 
this problem. 

References 

[1] Kurasov V.B., Physica A 226 (1996) 117 

[2] Kuni F.M., Grinin A. P., Kurasov V.B., Heterogeneous nucleation in va- 
por fiow. In: Mechanics of unhomogenenous systems, Ed. by G.Gadiyak, 
Novosibirsk, 1985, p. 86 (in Russian) 

[3] V.B. Kurasov Universality in kinetics of the first order phase transitions, 
SPb, 1997, 400 p. (in English) 

[4] Kurasov V. Phys. Rev. E, vol. 49, p. 3948 (1994) 

[5] Grinin A. P., P.M. Kuni, A.V. Karachencev, A.M.Sveshnikov Kolloidn. 
journ. (Russia) vol.62 N 1 (2000), p. 39-46 (in russian) 

[6] Grinin A. P., A.V. Karachencev, Ae. A. lafiasov Vestnik Sankt- 
Peterburgskogo universiteta (Scientific journal of St. Petersburg univer- 
sity) Series 4, 1998, issue 4 (N 25) p. 13-18 (in russian) 

[7] Grinin A. P., Kuni P.M., Sveshnikov A.M. KoU. Journ., 2001, volume 
63, N6, p.747-754 



26 



[8] Kurasov V.B., Deponed in VINITI Manuscript 2594B95 from 19.09.95, 
28p. 

[9] Kurasov V., Preprint cond-mat@xxx.lanl.gov get 0410774 
[10] Kurasov V., Preprint cond-mat2xxx.lanl.gov get 0410043 
[11] Kurasov V., Preprint cond-mat@xxx.lanl.gov get 0410616 
[12] Kurasov V., Preprint cond-mat@xxx.lanl.gov get 0412141 



27 



